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We construct semiclassical solutions of the symplectically covariant 
<**] ' Schrodinger phase-space equation rigorously studied in a previous paper; 

| we use for this purpose an adaptation of Littlejohn's nearby-orbit method. 

We take the opportunity to discuss in some detail the so fruitful notion of 
squeezed coherent state and the action of the metaplectic group on these 
states. 
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Abstract 



H ! 1 Introduction 



Let ^ be a square-integrable solution of Schrodinger's equation 

and set \l/ = U^tp, where the wave-packet transform U$ is defined as follows [5]: 
U^{x,p) = i^-Y' 2 e*v* [ e -**"*V(x')^(z - x')cTx' (1) 



here <p is an arbitrary rapidly decreasing function. The function '5 satisfies the 
equation 

ih— = H ph V (2) 
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where the operator i? p h is defined by the intertwining formula 

U$H ph = H(x, -ihd x )U,p; 
A straightforward calculation shows that we have 

Uj,(xip) = (\x + iHdp)U^ , U^{-ihd x ^) = (\p - ihd x )U<j>ip; (3) 
these relations motivate the notation 

-ffph = #(| x + ihd pi \pj - ihd x ) 
and we may thus rewrite J2J) as the phase-space Schrodinger equation 
d 

ih—^^H^x + ihdp^pj-ihdx)^. (4) 

In a recent paper [H] (also see jZj) one of us exposed in some detail the prop- 
erties of the solutions of this equation, originally proposed by Torres- Vega and 
Frederick |21l l2"2"j in the particular case where (f> is a Gaussian. 

In this paper we pursue the study of the properties of the solutions of this 
equation, and, in particular we show how to construct in an easy way approxi- 
mate solutions to that equation for a large class of initial wave functions. Our 
approach is an adaptation to phase space of Littlejohn's [17j nearby orbit method 
(we remark that Heller ^] has considered a more restrictive notion) . 



Notation 

The position vector will be denoted by x = [x\, x n ) and the momentum 
vector by p = (p±, ...,p n ), and we write z = (x,p) for the generic phase space 
variable. We will use the generalized gradients 



d - -2- 



3 - JL JL 

' U P — \ dpi ' d Pn 



and d z = (d x ,d p ). 

The symplectic product of z — {x,p), z' — (x',p') is denoted by a(z, z'); by 
definition: 

a(z, z 1 ) = p ■ x' — p' ■ x 
where the dot • is the usual (Euclidean) scalar product. In matrix notation: 

a(z,z) = (z ) Jz , J = 

y -*nxn Unxn 

The corresponding symplectic group is denoted by Sp(7i): the relation S G 
Sp(n) means that S is a real 2n x 2n matrix such that a{Sz, Sz') — a{z,z'); 
equivalents S T JS = SJS T = J. 

For zq = (xo,po) the Heisenberg-Weyl operator T(zq) acts on functions of 
x via 

f{z )ip(x) = e^P°- x -^ - X0 ^{x - x Q ). (5) 
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The Wigner-Moyal transform of a pair (ip, <fr) of functions in the Schwartz space 
«S(R£) is defined by 

WW, <j>){z) = {^f J e-i p - y 4>(x + iy)<p( x -y)d n y. (6) 

2 Coherent States in Phase Space 

In this Section we review known results about coherent states; the action of the 
metaplectic group on these states will be studied in the next section. 
We denote by <p h the standard coherent state: 

A*) = (^) n/ V*«' ; (7) 

for z — (xq,Pq) we define more generally the coherent state <p^ o = T{z )(f) h that 
is 

4> h zg { x ) =e*<-P°- x -?P°- Xa U h (x-x ). (8) 

(These states are often denoted respectively by |0) and \zo) in bra-ket notation 
|17|). The following properties are well-known (see for instance. [1511201 : for the 
sake of self-containedness we however give a proof. 

Proposition 1 (i) We have 

(a^s)"/ <t>l(x)$Md 2n z =5(x-y). (9) 
(ii) For every ip £ L 2 (R n ) 

and 

IMl£» = (a£s) n / l(V>,<M 2 d 2 "z . (11) 

(Hi) Let A be a continuous linear operator <S(R") > <S'(R™) with Schwartz 

kernel K; we have 

K ^y) = (^n) n J UU^)<M)d 2n ^ (12) 

2.1 Squeezed coherent states 

Despite the interest of the standard coherent states it is often advantageous to 
work with their generalization, called "squeezed coherent states" (they anyway 
appear automatically when one lets the metaplectic group act on coherent states, 
see Subsection 13.211 . A squeezed coherent state is a Gaussian of the type 

tfx,Y)(*) = (±) n/4 (^X) 1/4 ^ iX+lY)x2 (13) 
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where X and Y are real symmetric n x n matrices with X > 0; we have 
IIWx,y)IU a = We will find it convenient to set 

M = i(X + iY) , X = X T > , Y = Y T . 



and to write <p^ x Y \ = 4>m\ thus: 

rt 1 ( X ) = (±) n/i (detX)^e^* 2 , X = ImM. 
The Wigner transform of 4>m ^ s gi ven by the formula 

W^ M (z) = (^) n (detX)-^e-^ 2 (14) 
where G is the symmetric matrix 

an essential remark is that G is in addition symplectic (this fact was apparently 
first observed by Bastiaans |T|). More precisely, we have G — S T S with 

as results from a direct calculation; notice that S belongs to the isotropy sub- 
group of the Lagrangian plane f p = 0x W l in Sp(n). For z £ R 2n we set 



the Wigner transform of cf> Zo M is given by 

W4> h Z0 , M (z) = W<j> h M (z-z y, (17) 

it is thus a Gaussian centered at Zq. 

Let us generalize formula l|14f> by calculating the Wigner-Moyal transform 
^{4>zo,Mi 4>z' ,M'} °^ a P a ^ r °^ scmeeze( i coherent states; recall for this purpose 
the "Fresnel formula" 

(_i_)"/ 2 j e -^e~^ Kx ' 2 d n x' = (det K)-V*e-& K (18) 

valid for £ € C n , X = K T , Imif > 0; here (det X)- 1 / 2 = A~ 1/2 • • • A~ 1/2 where 
Xj 1//2 is the square root with positive real part of the eigenvalue XJ 1 of K^ 1 
(see HQE]). 

Proposition 2 We have 

W{^ M A h Ml ){z) = (±) n (detXX>)- 1/4 e-i Fz2 (19) 
where F is the matrix 

( 2iW(M-M>Y x M -i(M + W)(M -W)- 1 } 

\-i(M -Wy^M +W) 2i(M~W)- 1 )' ( ' 
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Proof. Setting Wm,M' = ^(Ai we nave 

W MM ,{z) = (^) n J + \v)~4k,(x- \v)<Ty 

that is, replacing M and by their expressions, 

Wm,m>{z) = C{X,X')e^ {M+JF)x2 J e-* py e^ x ' y) d n y 

where M = i(X + iF), M' = i(X' + iY'), the terms C(X, Y) and $(x, y) being 
given by 

C*(X,Y) = 2~ n (^-) 2n (detXX') 1/4 

y) = M (x + \y) 2 + W{x - \yf . 

Expanding &(x, y) and using (|18|) we finally get 

W MM ,{z) = (^{detXXr^e-r^ 

where F is the matrix ill'Ot (the left-upper block is obtained by using the identity 

M + W- (M - W)(M + M 7 )- 1 ^-/ — M 7 ) = 4M 7 (M + TF)~ 1 M. ) 

■ 

Noting the following formula, which is an easy consequence of the definitions 
of the Wigner-Moyal transform and the Heisenberg-Weyl operators: 

Lemma 3 For all zo, z' Q in M. 2 . 71 and f,g in L 2 (M. 2n ) we have 

W(T(z )i/j,T(z' )i>')(z) = e i(^o-z' ,z)-^( Zo ,4)) 

xW(ip,ip')(z-l(z +z' )) 
we immediately get the following generalization of Proposition [21 
Corollary 4 We have 

tn<,M,<,M<)(2) = (^)"e^ CT ^^^-^(-o,4))(detrf)- 1/4 e -* F(2 -5^- 
where F is given by \2U\) : hence in particular 

W^u^mX*) = (^)"(detXX')- 1 / 4 e-*^-- ) 2 . 

2.2 Phase-space coherent states 

For each <f> G the operator 

U^(z)= (f ) n/ V(V,0)(|z) 
is an isometry of L 2 (M") onto its range it follows that: 



Proposition 5 Let $JL = [/^(^ ). For eac/i ^ € TL4, we have 

*(*) = (ikT /(^*t)^tWfto (21) 

H*IU 2 = (2^)"/ K*,^)^! 2 ^* (22) 

Proof. Let ?/> be defined by '5 = C/^V; bi view of part (ii) of Proposition ^ 
(formula (|10|l ) we have 

hence, since U$ is continuous, 
In view of formula l)42|) we have 

(*,^ ) L2 - {u^u+tpjv = MO 

hence (12 1 B : formula (|22(1 follows by a similar argument from formula (|1 lfl . ■ 

3 Metaplectic Group and Coherent States 

The metaplectic group Mp(n) is a faithful unitary representation of Sp 2 (n), the 
double cover of the symplectic group Sp(n). There are several different ways to 
describe the elements of Mp(n) (see for instance 0^1113!); f° r our purposes 
the most adequate definition makes use the notion of generating function for 
a free symplectic matrix because it is the simplest way to arrive at the Weyl 
symbol of metaplectic operators (and hence to their extension to phase space). 

The interest of the metaplectic representation comes from the fact that it 
links in a crucial way classical (Hamiltonian) mechanics to quantum mechanics 
(see for instance ^S] or [5| and the references therein). Assume in fact that H 
is a Hamiltonian function which is a quadratic polynomial in the position and 
momentum variables (with possibly time-dependent coefficients): thus 

H( Z> t) = ±z T H"(t)z=±H"(t)z 2 

where H"(t) is a symmetric matrix (it is the Hessian matrix of H). The associ- 
ated Hamilton equations z = d z H(z,t) determine a (generally time-dependent) 
flow consisting of symplectic matrices St- We thus have a continuous path 
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t i — ► St in the symplectic group Sp(n) passing through the identity at time 
t = 0: Sq = I. Following general principles, this path can be lifted (in a unique 
way) to a path t i — > St in Mp(n) such that So = I (the identity in Mp(n)). 
Choose now an initial wavefunction in, say, <S(R") and set ij)(x, t) = Stip (x). 
The function ip satisfies Schrodinger's equation 

d'ip 

ift— = H(x 1 ihd x ,t)ip 
where the operator H(x,ihd x ) is given by 

H{x,ihd x ,t) = ^(x,ihB x ) T H"{t)(x,ihd x ). 

3.1 Description of Mp(n) and IMp(n) 

Let W be quadratic form of the type 

W(x, x') = \Px 2 - Lx ■ x' + \Qx 2 

with P = P T , Q = Q , det L ^ 0; (we have set Px 2 = Px ■ x, etc.). To each 
such quadratic form we associate the generalized Fourier transform 

S w ,m^(x) = (4m) n/2 i m \fttetL\ J ei w ^ x '^(x')d n x' (23) 

where m corresponds to a choice of the argument of det L modulo 2tt. One 
proves (see |3]) that every S € Mp(n) can be written (of course in a non- unique 
way) as the product of two operators of the type l|23ll , and that the integer 

m[S) = m + m' - Inert (P' + Q) (24) 

is independent modulo 4 of the factorization Sw,mSw',m> of S; the class modulo 
4 of m(S) is the Maslov index of the metaplectic operator S. Since Mp(n) is a 
realization of the double cover of Sp(n) there exists a natural projection 

tt Mp : Mp(ra) — > Sp(n); 

that projection is a 2-to-l group epimorphism defined by the condition: Sw = 
7r Mp (Sw,m) ^ the free symplectic matrix generated by W, that is: 

/ \ a r r r\ , . / P = d x W(x,x'), 

(x,p) = S w (x >P )^i p > = - dxl w{x,x>). 



Writing S € Sp(n) in block form: 

s ={ct) < 25 > 

S is free if and only if det B ^ 0, and the generating function W is the quadratic 
form in 2n variables 

W(x, x') = \DB- X x 2 - {B- X fx ■ x' + \B- x Ax 2 (26) 

We will use the following covariance properties of the metaplectic operators: 
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• We have 

ST(z ) = T(Sz )S (27) 
for all S £ Mp(n) and z £ R 2n ; 

• If S e Mp(ra) has projection S = tt Mp (S) then 

W(S^,S(f>) = W(4>,(t))oS- 1 . 

The group generated by the Weyl-Heisenberg operators T(zq) and the op- 
erators S £ Mp(n) is a group of unitary operators in L 2 (W l ); it is denoted by 
IMp(n) and called the inhomogeneous metaplectic group. 

3.2 The action of Mp(n) on coherent states 

To describe the action of S £ Mp(n) on the squeezed coherent states <j) Zo M we 
will need the following Lemma. Let us denote by £(n) the Siegel half-space, 
that is 

E(n) = {M : M = M T , ImM > 0} 
(M is a complex n x n matrix). 

Lemma 6 Let S £ Sp(n) 6e #wen 6y andM e £(n). J7ien det(A+BM) ^ 
0, det{C + DM) ^0 and 

a(S)M =(C + DM) (A + BM)~ l £ £(n) (28) 

(in particular a(S)M is symmetric) , and 

a(SS')M = a(5)a(5')Af. (29) 

XTie action Sp(n) x E(n) — ► £(") defined by is transitive: if M, M' £ £(n) 
then there exists S £ Sp(n) sucft £/ia£ A/' = a(S)M . 

For a proof see 0]^]; also see the preprint by Combescure and Robert. 
Notice that (|29fl implies that S i — ► en (5) is a true representation of the sym- 
plectic group in the Siegel half-space. 

Let S £ Mp(n) have projection 



S = TT MP {S) = 



A B 
C D 



on Sp(n); then 



1 -m(S) 

5^ ft (x) = 4 ; ex P 

1 ttW v /ditp + iS5 



where the branch cut of the square root of det(A + iB) is taken to lie just under 
the positive real axis; m(S) is the Maslov index l|24(l of S. Using this formula 
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S4Q {x) is easily calculated: since by definition = T(zo)4> h the metaplectic 
covariance formula 127|) immediately yields 

s 4>z { x ) = T(Sz Q )S(f> h (x) 

(see Littlejohn [T2| for an explicit formula). 

The results above can be generalized to arbitrary squeezed coherent states: 

Proposition 7 Let <f> M , M G be a squeezed coherent state and S G 

Mp(n), S* = vr M P(5). We have 

S&M = <t>^(S)M , S<t>z ,M = T (Sz())(j)l { s)M- (30) 

(see for instance pEfllT7|). 

The Gaussian character of a wavepacket is thus preserved by metaplectic op- 
erators; as a consequence the solution of a Schrodinger equation with Gaussian 
initial value remains Gaussian when the Hamiltonian operator is associated to 
a quadratic Hamiltonian function. 

3.3 The Weyl symbol of a metaplectic operator 

For S S Sp(n) such that det(5 - I) ^ we set 

Ms^yis + rKs-i)- 1 

(it is the symplectic Cayley transform of S) . In one of us proved the following 
result: 

Proposition 8 (i) If S\y, m * s such that Sw — ^ Mp (Sw,m) has no eigenvalue 
equal to one, then 

- / 1 \ " Inert W xx r 

**"*<*> = (ss) y fl*«s-»\ J'*'" - ° n '' W * )i '"* (31) 

for ip G S*(i?™); here Inert W xx is the number of negative eigenvalues of the 
Hessian matrix of the function x i — > W(x,x). (ii) Every S G Mp(n) can 
be written as the product of two metaplectic operators of this type. (Hi) More 
generally, every S G Mp(n) with projection S such that det(S — I) ^ can be 
written as 

/ i \« ■!/(§) r 

S * M = (ss) 7mW=Wi I ^ M "" T ^^ < 32 > 

where v(S) is the Conley-Zehnder index of S. 
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Formulae 1|31|) and (|32[1 yield the Weyl representations of the metaplectic 
operators; for a definition and a precise study of the Conley-Zehnder index 
v(S) appearing in see [§]• In [S] it was also proven that formula can 
be rewritten alternatively as 

^ = (i) il/(5 V|dct(5-/)| J e-i< Sz ^T{{S - I)z)d 2n z (33) 



or 



^ = (i) *' y(§ V|det(S-/)| J f(Sz)f(-z)d 2n z. (34) 

4 The Nearby Orbit Method 

We now have the background material that is needed to study the nearby-orbit 
method for the phase-space Schrodinger equation. We begin by describing this 
method in the usual case; we refer to Littlejohn's excellent review ^7] for a 
discussion and interpretation. 



4.1 The "standard" case 

Consider an arbitrary (possibly time-dependent) Hamiltonian function H on 
R^ n ; we denote by (ft) its flow: t i — ► ft(zo) is the solution of Hamilton's 
equations z = Jd z H(z,t) passing through the phase-space point zq at time 
t = 0. We have, for every z, 

ft(z) = ft(z Q ) + S t (z )(z - zq) + 0((z - z ) 2 ) 

where St(zo) = Df t (zo) is the Jacobian matrix of f t at zq. Since Hamilto- 
nian flows consists of canonical transformations, St(zo) is a symplectic matrix: 
St(z~o) S Sp(n). Suppose now the point z is close to zq] then 

ft(z) « ft(zo) + S t (z )(z - zq); 

the nearby orbit approximation consists in replacing ft{z) by ft{zo) + St(zo)(z — 
z ). Noticing that 

ft(z ) + S t (z Q )(z - z Q ) = T(f t (zQ))St{zQ)T(zQ)~ 1 Z 

where T(zq) is the translation operator z i — > z + zq, an educated guess is that we 
will obtain the quantum counterpart of this approximation by replacing trans- 
lations by Weyl-Heisenberg operators, and symplectic matrices by metaplectic 
operators. In fact this guess is correct up to a phase factor arising from the 
non-commutativity of the Heisenberg-Weyl operators: the "true" semiclassical 
approximation to the solution ip of Schrodinger's equation 

ifi^-=H^ , V^ = 0) = Vo (35) 
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is given by 



U(t,zo)ip = e" 



^-'>T(/ t (zo))£ t (zo)T(z rV, 



o 



(36) 



where the phase ^(t) is given by 



70o,i)= f ±a{z(t'),z(t')) - H(z(t'),t')dt' 



(37) 



o 



One shows (j3J El) that the function ip — U(t,zo)ip is the solution of the 
Schrodinger equation associated to the Hamiltonian function 



obtained by truncating the Taylor series for H at ft(zo) after second order. 
The approximation (|36|) is valid when the initial state ip is localized near 
zo = (xq,po), that is when when xq and po are chosen to be the position 
and momentum expectation vectors (xiP ,iPq)l 2 and po = (pip ,ip )L 2 a t initial 
time t = 0; the metaplectic operators St(zo) ar e obtained by the usual lifting- 
procedure: the function t 1 — > St(zo) is a path in the symplectic group passing 
through the identity of Sp(n) at time t — 0; the function t i — > St(-Zo) is then 
the path in the metaplectic group passing through the identity at time t = 0, 
and such that ir Mp (St(zo)) = <St(,Zo) for all t. Formula l|3t)|) has the following 
interpretation (Littlejohn ^7]): the operator T(zq)^ 1 rigidly translates the ini- 
tial wavepacket to a wavepacket centered at the origin; the metaplectic operator 
St(zo) then "squeezes" it, and T(/ t (zo)) finally rigidly translates this squeezed 
wavepacket so that it becomes centered at ft(zo)', one then multipliers the result 
by e^ 20 '*). 

Let us work out explicitly formula l|36l) when the initial wavepacket ip is a 
squeezed coherent state: 

Proposition 9 Let <j)^ o M be an arbitrary squeezed coherent state. We have: 



H Zo (z,t) = H(f t (z ),t) 



+ H'(f t (z ),t)(z - f t (zo)) + -H"(f t (z ),t)(z - f t (z )) 2 



U(t,z )cl>l M = ei^ t Ui (zo)M s t (z ))M 



where a(St(zo)) is defined by $28\) with S = St(zo). 
Proof. By definition 4> h ZQ M = T(zo)0m hence 

U(t,z )4> h ZoM (x) = e***'*>T(/ t (^))S t (*b)^t 



that is, taking the first formula (|30[l into account, 



t/(Mo)<M = e^ z ^fUt{z ))€ {St(zQ )) 



M 
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which completes the proof. ■ 

Formula (UUil) can be extended to arbitrary square-integrable initial wave- 
functions using the machinery of coherent states reviewed in Section[21 applying 
formula JTDJ to i(j £ L 2 (R") one has 

and one then takes as semiclassical solution 

u(t)M*) = (^)"|w„,<)L^(Mo)<Wrf 2n ^. 

Since by definition T(zo)^ 1 (f>^ = r 1S the standard coherent state, this formula 
can be rewritten 

UMoW = (271)"/ (^ I <) i2 T(/ t U))e^ ( * ) 5 t (zo))0 ;i ( a :)d 2 ^o (38) 
(cf. equation (7.33) in fT7]1. 

Remark 10 Observe that neither U(t,Zo) nor U(t) are linear operators. 

Of course, a natural question is arising at this point, namely "How good are 
the semiclassical approximations (j2HJ, IpTfjl ?". Very precise estimates have been 
given in [D] ( see the discussion in the last section) ; we will content us 
here with quoting the following result: 

Proposition 11 Letifi be the exact solution to Schrddinger's equation J^35)) with 
initial datum ip = (j) ZQ . Then for each T > there exists Ct > such that 

\\U(t,z )cf> h Zo - i/)(;t)\\v<K) ^ (39) 
provided that the Hamiltonian function H satisfies uniform estimates of the type 

\d?H(z,t)\<C' atT (l + \z\r (40) 
for \a\ > m, \t\ < T and z e K 2n . 

(It is in fact possible to obtain precise bounds for the constant Ct; see [H])- 
4.2 Weyl operators on phase space 

In [H] (also see [0] for details) we noticed that the wave-packet transform Q is 
related to the Wigner-Moyal transform @ by the simple formula 

U^{z) = (f ) n/ V(V>,0)(iz) , ^ei 2 K) (41) 

It follows that we have 

(U^t/j, [/^')£ 2 (Kf* ) = (^VOl^rs) (42) 
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and hence each of the linear mappings is an isometry of L 2 (R") onto a closed 
subspace H$ of L 2 (R 2 ™) (the square integrable functions on phase space). It 
follows that UiUtf, is the identity operator on L 2 (R™) and that = U^U^ is 

the orthogonal projection onto the Hilbert space Ti.^,. Defining T p h(^o) by 

T ph (z )*(z) = e-TM^)m( z _ ZQ ) (43) 
we have the fundamental relation 

f ph (z o )f/ = Uffizo). (44) 
Formula l|44|l allows us to associate to every Weyl operator 

&K X ) = (2^)"/ a a {z a )T{z a )ij{x)d 2n z 
the phase space operator 

4**0*) = (2^)" j ' a a {z )T pYl (z^{z)d 2n z ; (45) 

of course 

A ph C/ = U+A. (46) 

Formulae Ij45(l and H31(l allow us to give an explicit description of the action of 
the metaplectic representation on phase-space functions: if det(S' - I) ^ wc 
define 

/ 1 \ n/2 j"( s ) r 2- 

s ^>=(^) vmf^ J'~""' TMz ° mz)d2nzK <47) 

the operators 5 p h are in one-to-one correspondence with the metaplectic oper- 
ators S and thus generate a group which we denote by Mp p h(n); that group is 
of course isomorphic to Mp(n). Of course, as an immediate consequence of The 
following equivalent formulae are immediate consequences of the corresponding 
expressions (|3*3*| and (|3^|l of S: 

•V = (£)V (s V|det(S-J)| J e-^ Sz ^f ph ((S~I)z)d 2n z (48) 

and 

S P h = (^)'V (S V|det(S-/)| J f ph (Sz)f ph (-z)d 2n z. (49) 

Notice that the well-known "metaplectic covariance" relation A o S = S'^^AS 
valid for any S € Mp(n) with projection S € Sp(n) extends to the phase-space 
Weyl operators Ap^: we have 

SphT p h(z )S~^ =f ph (Sz) , A o S p h = S^AphSph- (50) 

In Subscction l2.2l we defined coherent states in phase space. The metaplectic 
action on coherent states described in Proposition carries over to this case 
without difficulty, yielding the formulae 
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4.3 Nearby-orbit method in phase space 

Let us now state and prove the main result of this paper: 

Proposition 12 (i) LetU(t, zq) be the semiclassical propagator for Schrddinger's 
equation and set ip — U(t, Zo)ip . The wavepacket transform U<f, takes ip to the 
function defined by 

* = e^^hiMzoMStizo^hfizo^o (51) 
with = U^ipQ. (ii) In particular, if ip = (f> 

U^(U(t, z )<) = e^ z ^f ph (f t (z ))(S t (z )) ph ^. (52) 
Proof. Set ip — U(t, zo)ip ; by definition of U(t, zq) we have 

$ = e *7(*o,t)f (f t ( Zo ))S t (zo)T (zoT^o 
hence, using successively and l |3o) . 

= e^( Zo ^C/4f(/ t (z ))^(2o)r(2o)-Vo 

= e ^^' t )f ph (/ t (z ))(^(z )) ph f(z )p h 1 ;7^ . 

Formula H52[l follows since we have 

nz )^ h Z0 = Tiz^uMo) = Umzo)- 1 ^) = U^ h ) 

and = U${<t> h ). ■ 

The result above therefore suggests that the phase-space version of the semi- 
classical nearby-orbit propagator U(t, zq) should be given by the formula 

U ph (t, z ) = e^ (z ^f P Mt(zv))(S t (z )) ph f (z ) ph 1 (53) 

and we have: 

Proposition 13 Let '5 be the solution of the phase-space Schrodinger equation 

a* 

iH—=H ph ^ , tf(t = 0)=$* 

with = U^,((p z ). Suppose that H satisfies the conditions in Proposition 
Then, for \t\ < T there exists a constant Ct > such that 

\\U ph (t,z )$ h Zo ~ *(.,t)|| L2(Rr) < C T Vh. (54) 

Proof. The inequality Q54J) is an easy consequence of the estimate (j2HJ): the 
solution iff is given by — U<f,(ip(.,t)) where ij; is the solution of the usual 

Schrodinger equation 

iH^=Hi>, ^(i = 0) = <; 
since is a linear isometry we have 

||^ph(t,abX-*(-,*)||^cRS-) = II^Mo)^-^-,*)!!^) <C T Vh. 



14 



5 Conclusion and Discussion 



There are several problems and questions we have not discussed in this paper, 
and to which we will come back in forthcoming publications. Needless to say, 
there is one outstanding omission: we haven't analyzed the domain of validity of 
the nearby-orbit method much in detail. There are many results in the literature 
for the standard nearby-orbit method. For instance, Hagedorn ^3 ^] obtains 
precise estimates using the Lie-Trotter formula; similar results were rediscovered 
and sharpened by Combescure-Robert [3] using the Duhamel principle. Also see 
[TU] (Ch.2, §2.1). It shouldn't be too difficult to obtain corresponding estimate 
for the phase-space Schrodinger equation, using the properties of the wavepacket 
transform. It is on the other hand well-known that there are problems with 
long times when the associated classical systems exhibits a chaotic behavior; as 
Littlejohn points out in |17j . the nearby orbit methods probably fails for long 
times near classically unstable points. 
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